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Polariton lasing in a multilevel quantum dot 
strongly coupled to a single photon mode 
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We present an approximate analytic expression for the photoluminescence spectral function of a 
model polariton system, which describes a quantum dot, with a finite number of fermionic levels, 
strongly interacting with the lowest photon mode of a pillar microcavity. Energy eigenvalues and 
wavefunctions of the electron-hole-photon system are obtained by numerically diagonalizing the 
Hamiltonian. Pumping and photon losses through the cavity mirrors are described with a master 
equation, which is solved in order to determine the stationary density matrix. The photon first-order 
correlation function, from which the spectral function is found, is computed with the help of the 
Quantum Regression Theorem. The spectral function qualitatively describes the polariton lasing 
regime in the model, corresponding to pumping rates two orders of magnitude lower than those 
needed for ordinary (photon) lasing. The second-order coherence functions for the photon and the 
electron-hole subsystems are computed as functions of the pumping rate. 

PACS numbers: 71.36. +c,42. 55. Sa,42.55.Ah 
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I. INTRODUCTION 

Excitonic polaritons are quasiparticles made up f rom 
strongly coupled electron-hole pairs and photon^i^. 
They are experimentally realized in semiconductor op- 
tical microcavities with embedded quantum wells. The 
small volume of the microcavity, high reflectivity of 
its walls, and quasircsonancc condition between the 
confined-photon and excitonic energies guarantee the 
strong coupling regime. 

At very-low excitation rates, in mean only a sin- 
gle quasiparticle lives inside the cavity. With increas- 
ing excitation power, however, an abrupt increase of 
ground-state occupation takes place due to the quasi- 
bosonic statistics of the polaritons. A threshold be- 
haviour of the photoluminescence is observed. This be- 
haviour has been int erpreted as Bose-Einstein condensa- 
tion ofpolaritons^El or as a dynamical effect (polariton 
lasing)^. The latter position is motivated by the exper- 
imental demonstration that thermalization mechanisms 
are not effective. 

In the present paper, we start from the idea of the po- 
lariton lasepSl, where pumping provides a reservoir from 
which the low-lying polariton states are populated. Un- 
like common lasers, no population inversion is required 
and the active medium (the cxcitons) is strongly inter- 
acting with the cavity photons, forming the quasibosonic 
polaritons. 

The theoretical description of polaritons faces the dif- 
ficulties inherent to a many-particle strongly-interacting 
system working under a non-equilibrium pumping 
regime. Our strategy to tackle this problem is based 
upon two simplifications. First, we consider a finite 
systerrPEla that is a single photon mode, and a fi- 
nite number of single-particle states (ten) for electrons 
and holes. Then, the electron-hole-photon many-particle 
Hamiltonian is numerically diagonalized in order to find 



the energies and wavefunctions of the system. We stress 
that both Coulomb and electron-hole-photon interactions 
are treated exactly in our scheme. Second, we compute 
the stationary density matrix from a master equation 
which accounts for photon losses through the cavity mir- 
rors and pumping. The master equation is solved in a 
truncated set of many-particle states. Notice that these 
simplifications preserve the main ingredients of the prob- 
lem: the existence of fermionic and bosonic degrees of 
freedom, the strong coupling between them, the existence 
of a finite number of single-particle states for fermions 
(around 10^ in Ref. j5], 10 in our model) participating in 
the conformation of polaritons, a stationary state reached 
when pumping and losses are equilibrated, etc. 

Strictly speaking, our model describes a quantum dot 
supporting a few excitonic states and strongly interact- 
ing with the lowest photon mode of a thin micropillar. 
It covers a n intermediate region between the two-level 
dolj io | ii|i 2l and the infinite system (well) ^^\ It is simple 
enough to allow exact diagonalization but, at the same 
time, complex enough to capture many of the properties 
of the infinite system. 

The plan of the paper is as follows. In Sec. |TT] the 
model is described in details. In the next section, we 
briefiy sketch the algorithm for the numerical diagonal- 
ization of the Hamiltonian, and show a few results for 
the energy spectrum and matrix elements of operators. 



In Sec. IV we present the master equation for the density 
matrix and show typical occupations of many-polariton 
levels for low, intermediate, and relatively strong pump- 
ing rates. In Sec. |V] the way of obtaining the exact 
photoluminescence (PL) spectral function, and the ap- 
proximations leading to the simplified expression used in 
the paper are clarified. From this expression, we com- 
pute the intensity, position and linewidth of the main 
PL peak as functions of the pumping rate. Sec. |VI| is 
devoted to second-order coherence functions. Finally, in 
the last section, we summarize the main results of the 



paper. 



II. THE MODEL POLARITON SYSTEM 

As mentioned in the preceding section, we study a fi- 
nite polariton system. A GaAs micropillar with radius of 
about one micron or lower is considered, in such a way 
that the lowest photon mode is well separated from the 
higher modes , and we can assume that a single photon 
mode is coupled to the lowest electron-hole states. The 
active medium inside the cavity is described by a finite 
number of harmonic-oscillator states for electrons and 
holes, as shown in Fig. IT] The number of single-particle 
states (ten) is dictated only by practical reasons: the 
dimension of the many-particle Hilbert space grows ex- 
ponentially with the number of states. This finite system 
could be a good model for a quantum dot inside a thin 
micropillar, and even could be used to obtain the quali- 
tative behaviour of quantum well-based micropillars. 

The interaction Hamiltonian includes electron- 
electron, hole-hole, and electron-hole Coulomb interac- 
tions as well as electron-hole-photon coupling, the latter 
in the rotating- wave approximatiorP^ 
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The effective band gap, Egap, is taken as 1500 meV 
for GaAs. A is the detuning of the photon mode with 
respect to Egap. The harmonic oscillator energies are 
much smaller than Egap. We will neglect them in the 
single-particle energies of electrons and holes, and will 
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write: T> ' = Q, T> ' = Egap- g is the electron-hole- 
photon coupling strength. Notice that we are including 
only spin-up electrons, spin-down holes and one "circu- 
lar" polarization of photons in Eq. (fTj) . A model with the 
two photon polarizations, which, however, would dramat- 
ically increase the dimension of the Hilbert space, would 
make possible the study of interesting features such as the 
spontaneous build-up of coherence between "left-handed" 
and "right-handed" polaritongi^. /3 is the strength of 
Coulomb interactions, and (i,j||r, s) - the dimensionless 
matrix elements among harmonic oscillator states. 

The oscillator states are labeled by two quantum num- 
bers: the number of zeroes in the radial wave function, k, 
and the angular momentum projection along the cavity 
axis, I. The hole state, i, in Eq. M is the conjugate 
of the electron state i, that is, has the same k, but the 
momentum is — Z. This means that the photon interacts 




FIG. 1: (Color online) Schematic representation of the har- 
monic oscillator states for electrons and holes in the model. 



only with electron-hole pairs with zero angular momen- 
tum. As a consequence, the total angular momentum of 
the electron-hole system: 



^-EOr'+^ 



i(e) , ,(h) 



(2) 



is a conserved magnitude. In addition, the Hamiltonian, 
Eq. (fTj), preserves the polariton number: 



Npoi = N, 



pairs 



Nph = 



'e 



(4e^ + h\h^ + a) a. (3) 



We notice the similarity between ours and a finite 
Dicke model^^. The infinite Dicke model has been used 
to describe polaritons in microcavitiea^. The main dif- 
ference with our approach is the following. In the Dicke 
model of polaritons, we first solve for the excitons and re- 
tain only the ground state. Multiexcitonic states are not 
considered. This, may be, is a good approximation for 
far-appart, small (not supporting multiexcitons) quan- 
tum dots in a microcavity. 

Many-particle states with fixed Npoi and L are con- 
structed in the next section. We give here a preview 
in order to compare with the traditional picture of non- 
interacting polaritons. We take for the parameters the 
values, g — 'i meV, and (3 — 2 meV. The latter is a rea- 
sonable value for GaAs, leading to an exciton binding en- 
ergy of a few me Vs. The high value oi g is, however, not 
intended to be realistic. It is chosen in order to illustrate 
the interesting regime, not studied so far, where photon- 
pair coupling and Coulomb interactions are comparable. 
In Fig [2] we show all of the states with Npoi = 1 in the 
model. We joined with a dashed line the lowest states 
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FIG. 2: Polariton states with Npoi = 1. 



in each L tower (the yrast states) in order to conform 
the lower polariton (LP) branch. The upper polariton 
(UP) states, on the other hand, can be identified from 
the photoluminescence (PL) emission. We indicated in 
Fig. [2] the UP state in the L = sector. Notice that, 
because of the strong electron-hole-photon coupling con- 
stant, the UP state is pushed up to high energies in our 
model. In between LP and UP states there is a set of 
"dark" polariton states. They play an important role in 
the dynamics because they can not decay through pho- 
ton emission. Let us stress that, our high-5 regime could 
be of interest in other contexts, where ultra-high light- 
matter couplings have been reported^^. 

We shall see in Sec. |IV| that photon losses in the 
cavity and incoherent pumping can be modeled by two 
terms in the master equation for the density matrix. We 
will not include relaxation mechanisms inducing transi- 
tions between states in the same Npoi sector (acoustical 
phonons). As a result, the total angular momentum is 
conserved even when pumping and losses are taken into 
account. We will solve the dynamics in the L = tower, 
which will allow us to compute the PL emission along the 
pillar symmetry axis. 



Finally, let us comment about the truncation of the ba- 
sis of single-particle states in Fig. [T] For small quantum 
dots, this is a natural assumption. In thin micropillars, 
the number of states strongly coupled to the lowest pho- 
ton mode is large, but finite. In Ref. \5\ for example, 
it should be around 10**. In this sense, our model may 
be thought of as a scaled version of a micropillar. At 
larger excitation energies the electron-hole states behave 
incoherently and act as a reservoir for the lower polariton 
states. We partially take account of these higher excited 
states in our model of incoherent pumping (Sec. IV I. 



Coulomb interactions between polariton states and the 
reservoir, which is an additional source of decoherence, 
will be, however, neglected. 



III. EXACT DIAGONALIZATION RESULTS 
FOR THE ISOLATED SYSTEM 

For given Np^i and L, we diagonalize the Hamiltonian 
in a basis constructed from Slater determinants for elec- 
trons and holes and Fock states of photons. The wave 
functions are looked for as linear combinations: 



|/> 






Cs,,Sh,n\Se,Sh,n). 



(4) 



where Se and Sh are Slater determinants for electrons and 
holes with the same number of particles, Npairs, and the 
number of photons is n = Npoi — Npairs ■ When Npoi = 
there is only one state, the vacuum. When Npoi = 1 
there are 17 states with L = 0. One of them is the state 
with one photon (no pairs), and the remaining 16 states 
correspond to matter excitations (no photons), that is, 
all possible combinations of one electron and one hole 
states with total angular momentum equal to zero. On 
the other hand, there are 256 states with Npoi = 2, 1746 
states with Npoi = 3, etc. As Npoi increases, the num- 
ber of eigenstates of H rises, reaching around 18000 for 
Npoi > 10. We use Lanczos algorithms^ to obtain the 
energies and wavefunctions of the lowest states in each 
sector. 

We give in Fig. |3] a schematic representation of the 
ground state wave functions with quantum number L — 



0, and polariton numbers Np 



1 (case (a)), Npoi = 6 



(case (b)), and Npoi — 600 (case (c)). The detuning pa- 
rameter is fixed to A = —3 meV. This value corresponds 
to quasi resonance. Indeed, in the Npoi = 6 case, the 
distribution is peaked around Npairs — 3, whereas in the 
large- A'po; limit it is peaked around Npairs = 5, that is 
the mean occupation of fermionic levels is near 1/2. No- 
tice that the mean number of photons is around 595 in 
the latter case. 

In Fig. H (a) the many-particle effects on polariton 
(photon) emission are made evident. We plotted the en- 
ergy difference Egs{Npoi) - Egs{Npoi - 1) - Egap as a 
function of Npoi ■ A persistent blueshift towards the pho- 
ton energy (equal to A) is noticed as Npoi is increased. 

On the other hand, in Fig. Ill (b) the energy difference 
Eijp — Elp is plotted as a function of Npoi- For large 
Npoi numbers, this difference behaves like 2^ g^ Npoi- 

The obtained wave functions may be used to com- 
pute matrix elements of operators. As it will be seen 
in the next section, the most important matrix ele- 
ments related to photon emission and losses are {F\a\T)^ 
where the many-polariton states F and / are such that 
Npoi(F) = Npoi(I) - 1. We show in Fig. [5] (a) the matrix 
elements squared |(F|a|/)p for transitions from Npoi = 2 
states to the one-polariton ground state. A Lorentzian 
with r = 0.1 meV is used to smear out the transitions. 
The analogs of UP and LP states are also clearly dis- 
tinguished here and in any Npoi sector. The transfer of 
population from the UP state with Npoi polaritons to the 
LP state with Npoi ~ 1 polaritons will be a key ingredient 
in the dynamics, as will become clear in the next section. 
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FIG. 4: (a) Energy shift of photon emission from the ground 
state as a function of Npoi- That is, AE — Egs{Npoi) — 
Egs{Npai - 1) - Egap- (b) Scahng in the Eup - Elp energy 
difference. 
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FIG. 3: Weights of the sectors with given Npairs in the ground 
state wave functions. Case (a) corresponds to L = and 
Npoi = 1, case (b) to Npoi — 6, and case (c) to Npoi = 600. 



In Fig. 5] (b) we draw the absolute value of the matrix 
elements '(i^|a|/)| in the low- A'po; sectors. When Npoi > 
1, only the lowest 20 states are used to construct the 
matrix. Notice that the analogs of LP and UP states are 
always included among these 20 states. We computed the 
matrix elements for Npoi < 600 (a matrix of dimension 
around 12000) and stored them in a file. A second file 
contains the energy eigenvalues. They are the input files 
for the dynamics, discussed in the next section. 



IV. MASTER EQUATION DESCRIPTION OF 
PUMPING AND LOSSES 



The actual polariton system is not isolated. Photons 
escape mainly through the cavity mirrors. The spon- 
taneous pair decay through leaky modes of the cavity is 
much less importantP^l, and will be neglected. In order to 
maintain a mean number of polaritons in the cavity, the 
system should be continuously pumped. As mentioned 
before, pumping comes from excited pair states decou- 
pled from the photon field, which may decay through 



emission of optical phonons, for example. We will, how- 
ever, neglect the effects of Coulomb interactions with 
the excited pair states on the relaxation of the polariton 
states, and also neglect relaxation due to the emission of 
acoustical phonons by the polariton states, prevented by 
the bottleneck eSec^^, that is selection rules for energy 
and momentum of the transitions which can not be si- 
multaneously satisfied. These two effects, i.e. relaxation 
due to Coulomb interactions or to phonons, could be in- 
cluded in a latter stage, but in the present paper they 
will not be considered. This means that the density ma- 
trix of the polariton system should be determined from 

a dynamical equation. 

We will use a quantum dissipative master equatiorpSHH 
in order to describe photon losses and pumping: 



— = \H.p]-\ (2apa^ ^ a* ap — pa' a) 

at n ' 2 



P 



^(2cr|jpcr/j - aijaljp ~ paija\j), (5) 
i,j 

The parameter k accounts for photon losses through the 
cavity mirrors {hn w Egap/Q, where Q is the cavity qual- 
ity factor). In our calculations, we take k = 0.1 ps^^. 
Notice that k << g/h, thus our model system works 
under the strong light-matter coupling regime. On the 
other hand, the parameter P is a pumping rate. We will 
use a sort of homogeneous pumping, with equal proba- 
bilities for all states. To this end, we introduce lowering 
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FIG. 5: (Color online) (a) The matrix elements \{F\a\I)\^ for 
transitions from Npoi = 2 states to the one-polariton ground 
state. A Lorentzian with V = 0.1 meV is used to smear out the 
transitions, (b) The matrix |(F|a|/)| in the low-Npoi sectors. 
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and rising operators, crij\I) = \J), crjj|J) = |/), where 
Npoi{I) — Npoi{J) + l. As we are employing a finite num- 
ber of states (20) in each sector with given Npoi > 1, total 
pumping probabilities are finite. The absence of phonon 
thermalization is also the reason why L = states are 
decoupled from other states with L ^ 0. Thus, we will 
solve Eqs. ([5]) in the most relevant L — sector. In ad- 
dition, we will focus on the stationary solutions of Eqs. 
(l6|, that is, the l.h.s. of these equations equal to zero. 

The number of variables in Eqs. (l5| may be estimated 
as follows. In each sector with Npoi > 1 there are 20 
occupations, pu, and 20 x 19 = 380 coherences, ppi, 
with F ^ I. That is, 400 variables per sector. If we in- 
clude sectors with < Npoi < ^1"^^ ' the total number 



of variables is 400A'. 



pol 



When N^'^i'"'''^ = 10, for 



example, the system has 3998 equations. 

We solve the resulting linear system of equations for 
the stationary density matrix with A^p";"^ — 10, and 
found the remarkable fact that the coherences are three 
order of magnitude lower than the occupations, that 



FIG. 6: Occupations at three different pumping rates: low 
(upper panel), intermediate (polariton laser, central panel), 
and high pumping (lower panel). The detuning parameter is 
A = -3 meV. 



is the density matrix is approximately diagonal in the 
energy representation^. For example, for the set of 
parameters A = —3 meV, P = 0.01 ps^^, we get: 

E/^j|p/./l/E/P// = 7x10-4. 

In what follows, in order to extend the analysis up 
to relatively high polariton numbers {N^J""^ — 600), 
we will neglect the coherences. The number of variables 



reduces to 20iV^™°''^ 
stationary limit, Eqs. 



- 2. For the occupations in the 
Isl take the explicit form: 



= «5]|(/|a|J)r 



Pjj - npii 



Ek^h^)|' 



p 



E 



Pjj-PpnN^p{I). (6) 



N^al{J) = Np,i{I)~l 



where Nup{I) counts the number of states with polariton the auxUiary operator: 
number Npoi{I) + 1. We have 7V„p(l) = 17, Nup{I) = 20 



for 1 < Npoiil) < N^pTr^ and, finally, Nup{I) = for 

Npoi{i)^N^pTr^. 

The set of homogeneous linear equations (|6| should be 
complemented with the constraint. 



gaJJ={\J){I\{t + T)a{t)), 



(10) 



satisfies with respect to r the same master equation as 
the matrix elements p/j, with initial conditions: 



Y.pii^'^' 



(7) 



which corresponds to the conservation of probability. 

We show in Fig. [6] three regimes of pumping: low, 
intermediate, and large pumping rates, clearly differen- 
tiated by the patterns of occupations. In that figure, 
the y-axis corresponds to the occupations p//, whereas 
in the a;-axis the states are arranged in increasing order 
of the polariton number, Npoi- Recall that the first state 
is the vacuum with Np^i — 0, then we have 17 states 
with Npoi — 1, then 20 states with Npoi — 2, etc. The 
ground state in each sector with fixed Npoi is indicated 
by a square. 

At low pumping rates, the mean polariton number, 
defined as {Npoi) = Y.i PiiNpoi{I), is {Npoi) ~ 1- The 
state with the highest occupation is the vacumm. The 
ground-state occupations in sectors with Npoi > 1 are de- 
pressed. On the other hand, in the situation represented 
in the central panel of Fig. [6] {Npoi) is around four. The 
ground state occupations in sectors with Npoi < {Npoi) 
are enhanced with respect to the other states in the same 
sector. This is a kind of stimulated occupation of ground 
states. Finally, for large pumping rates the occupation 
in each sector with fixed Npoi is nearly uniform. In the 
example shown in the lower panel of Fig. m {Npoi) is 
around 24. A broad bell of occupied states ranging from 
Npoi ~ 12 to around 40 is observed. 

Once computed the stationary density matrix, one can 
estimate the photoluminescence response in the station- 
ary state. 



9a,ij\r=o = Xl^^l^^l^^ PKjit)- 
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In the stationary limit, t ^ oo, we get PKj{t) — Pjj^kj, 
and ga,ij{T -^ 0) = {I\a\J) Pjj- These initial condi- 
tions dictate that gajj behaves in the same way as the 
"vertical" coherences, that is, Npoi{I) = Npoi{J) — l. Re- 
call the equation for the vertical coherences, which may 
be obtained from Eq. (Is]): 



dr 



ga,ij = (iw/j - ^i.i)ga. 



IJ 



K,M 

I ^ {I\a^\M){M\a\K)go,Kj 

f E 9ajM{M\a^\K){K\a\J), (12) 



where uij = (Ej — Ei)/h, and 
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P . 



J2{\{K\a\I)\'+\{K\a\J)\'} 



K 



+ -^{Nup{I)+N^p{J)}. 



(13) 



The general solution of Eq. ( 12 1 is 



V. THE PHOTOLUMINESCENCE SPECTRAL 
FUNCTION 

In order to obtain the photoluminescence spectral 
function, S{ijj), we follow the lines sketched in paper [lO] . 
S{lo) is defined in terms of the first-order correlation func- 
tion of photons: 



S{uj) = -Re dTexp{~-iujT){a^{t + T)a{t)). (8) 



This function is to be computed with the help of the 
Quantum Regression TheorerrPS, which states that if we 
write: 



^{t + r)a{t))=Y,{J\a^\I)gojj, 



(9) 



9a,IJ 



Y.Cne''^^^u\ 



(14) 



where A„ and Xjj are, respectively, the eigenvalues and 
eigenvectors of matrix Bjj^mk, defined by the r.h.s. of 
Eq. ([12]), that is 



DlJMK^MK - ^n^IJ ■ 
M,K 



(15) 



I, J 



The coefficients Cn are determined from the initial con- 
ditions: 



J2CnX[-^^{I\a\J)p^-;\ (16) 

n 

The explicit expression for S{u;) is the following: 
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FIG. 7: (a) PL spectral function computed from the exact 
expression, Eq. (17 1, (b) The spectral function computed 
from the approximate expression, Eq. (191. The detuning is 
A = — 3 meV, and the pumping rate, P = 0.006 ps~^. 
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Where Djjn — (>/|a^|-f)C'„X|j , and the supraindexes r 
and i refer, respectively, to the real and imaginary parts 
of the magnitudes. The dimension of the matrix prob- 
lems given by Eqs. p^ and (fT6| is 17 + 20 x 17 + 20 x 



2). When TV^j^ = 10, for example, the 



20x(iv(7^ 
dimension is 3557. 

We notice that there is an approximate expression for 
S(oj) which is based on the fact that Ej — Ei ^ Egap ~ 
1500 meV (for GaAs), whereas Hk and HP are smaller 
than 1 meV. In a first approximation, we take only the 
diagonal terms in Eq. ( |12[ ), arriving to the following ex- 
pression for the correlation function: 



(a\t + T)a{t))\ 

^^|(/|a|J)|2 p(,f exp(*c.,,,-r,,,)r, (li 



I t 'rOO 

|2 ^(oo) 



/,,/ 



from which it follows that 



S{u) 



Y^ 



2 Joe)-, 



\{I\a\J)\'p^jyTjj 



(19) 
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FIG. 8: S{E) and the Lorentzian fit (dashed line) to the LP 
peak for the three cases of Fig. |6] 



The main advantage of expression ( 19 1 is the simplic 



ity. The luminescence from state J depends on the proba- 
bility, pjj, that the state is occupied, and on the matrix 
elements {I\a\J) for emission of a photon. The widths 
T jj have contributions from losses and pumping, the lat- 
ter is also a source of decoherence. 



The nondiagonal terms in Eq. (12) can only slightly 
modify the position of resonances, given by oj/j. They 
have a more appreciable effect on the widths. In Fig. 
n\ a comparison is made between the exact, Eq. (17 1, 
and approximate, Eq. (19 1, spectral functions. The pa- 



rameters are such that the mean number of polaritons 
is (Npoi) = 3.4. The lower energy emission has contri- 
butions from different peaks. We notice, by the way, 
that multimode emission in the polariton lasing regime, 
which is a manifestation of its non-equilibrium nature, 
has been nicely demonstrated recently^'^. The strongest 
peak, which we take as the definition of the LP, is more 
sharper in the exact scheme. We will, nevertheless, use 
expression ( 19 1 in order to obtain the behaviour of the 
PL even for very strong pumping rates {{Npd) « 500), 
where an effective weak coupling regime is established. 
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FIG. 9: Integrated intensity, position and linewidth of the 
lower polariton peak, coming from Eq. (19 1, as functions of 



the pumping rate. A = —3 meV. For the additional lines 
in the center and bottom panels see explanation in the main 
text. 



We fit the lower polariton peak to a Lorentzian (dashed 
line), from which the integrated intensity, peak position, 
and effective linewidth are extracted. In Fig. |8] we show 
S{E) and the corresponding Lorentzian fit for the three 
cases illustrated in Fig. [6] The main characteristics of 
the polariton emission are apparent in the figure. That 
is, a blueshift of the emission, and an increase of the 
linewidth as the pumping rate is increased. 

The upper panel of Fig. [9] shows the integrated inten- 
sity as a function of P for a detuning A — —3 meV. A 
threshold (change in the slope) at P« 3 x lO^-^ ps^^ is 
observed, corresponding to stimulated ground-state oc- 
cupations when the number of polaritons exceeds one. 
This is the "polariton laser" regime. At this threshold 
value, the peak position (center panel) begins a contin- 
uous blueshift towards the bare photon energy (1500- 





ji| 1 1 1 1 1 1 1 


1 


"1 1 '; 


^ , 


! 






w 


. 














■E 100 


r 




JK^'\ 


Z! 








JD 


! 






k 10 


- 




—M^ 










>. 


; 


m^ 


' 




. 


^*^ 


, 


1 1 




W 




r 


/ 


1 


03 




d 








^ 




c 


! 


f 


\ 


^ 0.1 


j»^ 


i 


-. 


c 


fil*'*'*'*^ 


1 




1504 


- 


1 




















■S 1500 


_ 




_ 


E 




r / 




'^ — - 


. 


/ f 


, 


>. 




/ r 




^ 1496 


- 


I / 
y 

t 
i/ 


- 


c 

LU 


-^ 


\ 


. 






i_ 






/ff 


0) 




1492 


__*^ 


1 .'mm 


- 




9:^^' 


1 1 


100 


P 


1 1^1 111 

o 


'V^, 




: 


"^ 


J^ 


> 


; 


nj 


jg^ 


03 


. 


o 


^^ 


^ 10 


- 


Q- 


^^ 


sz 


\ 


—J 


r^ 


-t-^ 




a a *^ 




"O 


. 


y^ 


. 


> 


■ 


/ 


■ 


^ 1 


_ 


{ 




r m^ 


p 


"1 


_l 


■ 1 


1 


■ ■1 1 



10"^ 10"^ 10"'' 10° 

Pumping (ps"^) 

FIG. 10: Same as Fig |9] for A = +5 meV. 



3=1497 meV), and the linewidth (bottom panel) starts 
increasing. In the "polariton laser" regime there is an in- 
terval where the linewidth saturates, and even decreases. 
This is better seen in the additional line (empty squares) , 
computed from the exact equations, Eqs. ( |T5p7 l. The 
decrease of the linewidth corresponds to maximum co- 
herence, as will become evident below. 

We draw an additional curve (dashed line) in Fig. [9] 
center panel, which refers to ground-state to ground-state 
transitions. This curve is constructed in the following 
way. For a given P, we find (Npoi). Then, the energy of 
the transition from the ground state of the system with 
polariton number equal to (Npoi) to the ground state of 
the system with polariton number equal to (Npoi) — 1 is 
found from Fig. ETa). Comparison with this curve shows 
that the excited states, and states with polariton number 
higher than the mean value determine the position of the 
LP peak. 

The right border of the polariton laser regime is con- 
ventionally set to P « 4 X I0~^ ps~^ in this figure. It is 





."1 
^il 1 


1 ' ' 

<Npol> 
<Npairs> 


1 1 ■; 


% 100 




CL 

z 

V 

% 10 




V 

1 






M| 1 


c 
o 

CL 


1 1 ' 


2.0 


9 ph • - 
-/21 n 


-^ 1.5 


g eh ° 
3lD DDDDDDDD 


1 .0 


III ■ 


----H □ B'"" 





10" 



10" 



10" 



lo"- 



Pumping (ps" 



FIG. 11: Upper panel: Mean number of polaritons and 
electron-hole pairs as a function of P. Lower panel: The 
second-order coherence functions at zero time delay for pho- 
tons and electron-hole pairs. The detuning parameter is 
A = -3 meV. 
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FIG. 12: Same as Fig. [ll]for A = +5 meV. 
rameters, but anyway they illustrate aspects of principle. 



VI. 



characterized by a second change of slope in the intensity 
curve, and a renewed increase of linewidth. Let us notice 
that, in this quasiresonant case, the mean occupation of 
fermionic levels becomes near one half in the border, a 
fact that could be appreciated below. For strong pump- 
ing rates, we observe a tendency to saturation in the 
position of the peak (towards the photon energy), which 
indicates the emergence of a new regime characterized by 
an effective, weak photon-matter interaction. 

In Fig. [To] we show results for positive detuning, 
A = +5 meV, where the electron-hole contents of the 
wave functions are higher. At this level, they look very 
similar to those reported in Fig. [9J We can approx- 
imately fix the limits of the polariton laser regime as 
9 X 10~^ < P < 6 X 10~^ ps~^. However, as it will be- 
come clear in the next section, the mean number of pairs 
is close to 5 near threshold, leading to a mean occupa- 
tion of fermionic levels, {Npairs)/^0, close to 1/2. This 
example shows that polariton lasing is not in antagonism 
with population inversion. Or, to set it in a different 
way, population inversion in these systems is not syn- 
onym of effectively weak pair-photon coupling. These 
results could be related to the small number of available 
states for fermions or the chosen values of the system pa- 



SECOND-ORDER COHERENCE 
FUNCTIONS 



In Figs. [TT]and 12 we show the mean number of polari- 



tons, (Npoi), the mean number of pairs, (Npairs), and the 
coherence properties of the photon and matter subsys- 
tems in the quasi-resonant and in the positive detuning 
case, respectively. We define the second-order coherence 
functions for photons and electron-hole pairs in terms of 
the one- and two-point correlation functions at zero time 
delay: 



(2) 

9ph 



'a*a* aa) 



(2) ^ {D^D^DD) 



(20) 



(21) 



where a is the photon annihilation operator, and D — 
J2i hjCi is the interband dipole operator. 

The coherence functions evolve from values larger than 
two at low pumping rates to values near one (Pois- 
son statistics, perfectly coherent state) immediately af- 
ter the threshold. Notice that the electron-hole sub- 
system reaches coherence more rapidly than photons 
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J2) 



J2) 



ideh "^ 9ph) possibly because of Coulomb interactions. 



(2) 



,(2) 



For large values of P, we get g),fl > g f^ ~ 1. Notice also 
that, in the resonant case, there is a pumping rate for 

(2) (2) 

which both g^^ and g^^ are approximately equal to one. 
This is the point of maximum coherence, and corresponds 
to a minimum of the linewidth. 

In the positive detuning case, the minimum of the 

(2) 

linewidth is reached at the point where 5^^ has a local 
maximum. For low pumping rates, (Npoi) and (Npairs) 
are very similar. They start differing precisely at the po- 
lariton lasing threshold, where the population of fermions 
is inverted. 



VII. CONCLUSIONS 

In conclusion, we have computed the stationary den- 
sity matrix, the photoluminescence spectral function, and 
the second-order coherence functions in a model polari- 
ton system describing a multi-level quantum dot strongly 
interacting with the lowest photon mode of a microcav- 
ity. The main features of polariton lasing, i.e. blueshift 
of the emission peak and increase of the linewidth as the 
pumping rate is rised, are reproduced by the model. Un- 
expected properties, such as the coexistence of polariton 
lasing and population inversion for positive detuning, are 
also manifested. 



Our polariton model with a finite number of degrees 
of freedom, could be positioned in between the two-level 
system, studied in Ref. [TOl and the infinite degrees of 
freedom systems considered, for example, in Refs. [T3] 
Our model is simple enough to be numerically diagonal- 
ized but, at the same time, complex enough to capture 
many of the features of the infinite system. In this sense, 
our results could be qualitatively compared with the ex- 
periment reported in Ref. |5J although the values of our 
model parameters are completely unrealistic. 

Finally, we should stress that we are aware of the lim- 
itations of our model. We understand, for example, that 
g^'^' does not rise further in the regime of strong pump- 
ing because the model does not include higher fermionic 
levels, which should become populated in this regime, 
that interact with the lowest polariton states, partially 
destroying coherence. 
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